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多 中 心 半 开 放 式 送 取 需 求 可 拆 分 的 车 辆 路 径 研究 


张 颖 钰 ， 吴 立 云 ” 
(河南 理工 大 学 工商 管理 学 院 , 河南 焦作 454003) 


摘 要 : 针对 多 中 心 半 开 放 式 送 取 需求 可 拆 分 的 车 辆 路 径 问 题 ， 构 建 了 以 车 辆 配送 距离 最 短 为 目标 的 多 中 心 半 开 放 
式 送 取 需 求 可 拆 分 的 数学 模型 。 设 计 大 变异 邻 域 遗 传 算法 进行 求解 ， 采 用 二 维 染 色 体 编码 及 顺序 交叉 策略 ， 同 时 
用 大 变异 策略 和 邻 域 搜索 策略 提高 算法 全 局 和 局 部 的 寻 优 能 力 ， 通 过 算 例 对 比 验证 了 所 提 模 型 与 算法 的 有 效 性 。 
例 实 验 表 明 ， 大 变异 邻 域 遗 传 算法 在 求解 多 中 心 物 流 配 送 车 辆 路 径 问 题 上 ， 求 解 质 量 较 优 、 求 解 效率 较 高 、 求 解 
果 较 为 稳定 。 同 时 还 验证 了 联合 配送 下 多 中 心 半 开 放 式 送 取 需求 可 拆 分 的 配送 模式 优 于 独立 配送 下 单 中 心 送 取 需求 
可 拆 分 的 配送 模式 。 研 究 成 果 不 仅 拓展 了 车 辆 路 径 问 题 ， 还 可 为 相关 快递 物流 企业 配送 优化 提供 决策 参考 。 
关键 词 : 车 辆 路 径 问 题 ; 多 中 心 ; 送 取 需 求 可 拆 分 ; 大 变异 遗传 算法 
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Multi-depot half open vehicle routing problem with split deliveries and pickups 


Zhang Yingyu, Wu Liyun’ 
(School of Business Administration, Henan Polytechnic University, Jiaozuo Henan 454003, China) 


Abstract: To solve the problem of multi-depot half-open vehicle routing problem with split deliveries and pickups, this article 
constructs a multi-depot half-open vehicle routing problem with split deliveries and pickups mathematical model with the 
shortest vehicle delivery distance as the goal. The solving process mainly use genetic algorithm based on big mutation and 
variable neighborhood search, which adopt two-dimensional chromosome coding and sequential crossover strategy, the big 


mutation strategy and neighborhood search strategy improve the global and local optimization capabilities of algorithm. The 
numerical examples verify the effectiveness of the model and the algorithm. Numerical examples show that the algorithm 
have benefits in solution quality, efficiency for solving the multi-center logistics distribution vehicle routing problem, and the 


result is relatively stable. Meanwhile, the multi-depot half-open vehicle routing problem with split deliveries and pickups is 
2 hs superior to the vehicle routing problem with split deliveries and pickups. The research results not only expand the vehicle 


routing problem, but also provide decision-making references for the distribution optimization of related logistics enterprises. 
re Key words: vehicle routing problem; multi-depot; split deliveries and pickups; genetic algorithm based on big mutation 


>< 0 引言 点 的 送 取 需 求 ， 有 助 于 提高 车 辆 载重 利用 率 ， 降低 物流 、 运 
GG 输 企 业 的 配送 成 本 。 目 前 ， 针 对 该 问题 的 研究 涉及 较 少 。 但 
四 近年 来 电 商 联盟 快速 发 展 ， 环 保 政策 相应 落实 ， 促 使 匡 其 包含 的 SVRPPD、MDHOVRP 等 问题 仍 是 当前 的 研究 热点 。 
= 家 逆向 物流 发 展 迅速 。 车 辆 路 径 问 题 (Vehicle Routing Problem,， 关于 SVRPPD、MDHOVRP 国内 外 学 者 都 对 其 模型 和 求 

VRP) 逐 渐 演化 成 联合 配送 下 多 中 心 带 回 程 的 运输 问题 , 其 中 解 方法 进行 了 广泛 和 深入 的 研究 。Mitra5] 首 次 深入 分 析 了 


包括 多 中 心 开 放 式 车 辆 路 径 问 题 (Multi-Depot Open VRP， SVRPPD 问题 ,并 构造 了 NH、SM 两 种 启发 式 算 法 对 其 求解 。 
MDOVRP)U、 多 中 心 半 开放 式 车 辆 路 径 问题 (Multi-Depot ”Archetti 等 人 MI 对 客户 点 需求 是 否 拆 分 作出 了 详细 的 科学 论 
Half Open VRP, MDHOVRP)D 以 及 多 中 心 闭 合式 车 辆 路 径 问 证 ， 实 验证 明 当 客 户 点 的 总 平均 需求 大 于 车 辆 容量 一 半 但 小 
题 。 随 着 客户 点 由 单一 的 配送 需求 演变 为 送 取 双 需求 ， 还 包 ”于 车 辆 容量 的 四 分 之 三 ， 并 且 需 求 方差 相对 较 小 时 ， 拆 分 可 
括 同 时 送 取 货 问题 (VRP with Simultaneous Delivery and 以 获得 最 大 的 收益 。 拆 分 带 来 的 收益 主要 体现 在 车 辆 配送 路 
Pickup, VRPSDP)61 和 送 取 需 求 可 拆 分 的 车 辆 路 径 问 题 (VRP ， 径 的 减少 ， 跟 客户 点 的 位 置 无 关 。Tang 等 人 中 在 解决 现实 中 
with Split Deliveries and Pickups, SVRPPD) 员 等 众多 的 拓展 问 第 三 方 物流 案例 中 , 采用 了 送 取 需求 可 拆 分 的 车 辆 路 径 模型 ， 
题 。 本 文 研究 的 联合 配送 下 多 中 心 半 开放 式 送 取 需 求 可 拆 分 ”并 加 入 了 时 间 窗 的 限制 , 运用 构造 型 启发 式 算法 -竞争 决策 算 
的 车 辆 路 径 问 题 (Multi-Depot Half Open VRP with Split ”法 对 模型 求解 .Nagy 等 人 [SI 采用 禁忌 搜索 算法 对 SVRPPD 进 
Deliveries and Pickups, MDHOSVRPPD) 是 综合 了 SVRPPD、 行 求解 ， 并 与 文献 [7] 的 算法 进行 了 对 比 ， 表 明 禁 忌 搜 索 算法 
MDHOVRP 及 联合 配送 特点 的 车 辆 路 人 径 问题 。 在 求解 精度 和 效率 上 均 更 优 。 对 于 MDHOVRP 的 研究 ， 

MDHOSVRPPD 问题 更 加 准确 的 描述 了 现实 物流 配送 的 复杂 ” ”Cordeau 等 人 四 设计 了 一 种 禁忌 搜索 算法 , 并 加 入 了 一 组 惩罚 
情况 。 如 在 一 个 多 配送 中 心 联合 配送 快递 的 区 域内 ， 配 送 车 ”函数 对 MDHOVRP 问题 求解 ，Ting 等 人 09 将 模拟 退火 算法 
辆 不 仅 要 配送 快递 ， 还 要 进行 快递 的 揽 件 操作 ， 当 客户 点 需 ”与 蚁 群 算法 相 结合 对 MDHOVRP 问题 进行 求解 ,就 勇 等 人 [1 
求全 部 满足 后 ,配送 车 辆 可 采用 就 近 原则 返回 任意 配送 中 心 ， 研究 了 带 时 间 窗 的 多 中 心 半 开放 式 车 辆 路 径 问 题 ， 并 设计 了 
共享 配送 网 络 。 半 开放 式 的 配送 中 心 及 合理 的 拆 分 末端 客户 ” 一 种 三 段 式 求解 方法 ， 第 一 阶段 先 用 K-mediods 方法 将 多 配 
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送 中 心 转 换 为 单 配 送 中 心路 径 , 第 二 阶段 用 GA-ACO 算法 进 Deliveries and Pickups, SVRPPD)。SDVRP 问题 可 以 看 做 是 
行 求解 ， 第 三 阶段 运用 节约 算法 将 单 配送 中 心 拆 分 成 多 配送 SVRPPD 问题 中 送 取 需求 其 中 一 个 为 0 时 的 特例 。 而 
中 心 ; 范 厚 明 等 人 (3 根据 蚁 群 算法 ， 改 进 信 息 素 更 新 方式 从 ”SVRPPD 问题 可 以 看 成 是 对 SDVRP 与 VRPSDP 的 扩展 ， 
而 提高 寻 优 能 力 , 对 MDHOVRP 问题 进行 求解 ; 王 勇 等 人 03] ”SVRPPD 的 示意 图 见 图 1(d)。 


将 遗传 算法 与 粒子 群 算法 结合 求解 MDHOVRP 问题 。 范 厚 明 5) 多 中 心 闭合 式 送 取 需 求 可 拆 分 的 车 辆 路 径 问 题 
等 人 [9 首次 研究 多 中 心 联盟 配送 与 同时 送 取 货 结 合 的 车 辆 随 着 现代 物流 业 的 发 展 ， 逐 渐 演 变 为 在 某 一 区 域内 有 多 
问 题 (Multi-Depot VRP with _ Simultaneous Deterministic 个 配送 中 心 ， 来 自 不 同 配送 中 心 的 车 辆 完成 所 有 客户 点 的 送 


货 任 务 后 返回 原配 送 中 心 ， 称 为 多 中 心 闭合 式 送 取 需 求 可 


Delivery and Stochastic Pickup based on Joint Distribution, 

MDVRPSDDSPJD)， 根据 问题 特征 ， 设 计 了 变 邻 域 文化 基因 拆 分 车 辆 路 径 问题 (Multi-Depot VRP with Split Deliveries and 

算法 求解 该 问题 。 Pickups, MDSVRPPD)。MDSVRPPD 问题 是 对 SVRPPD 问题 
综 上 所 述 , 现 有 文献 主要 集中 在 对 SVRPPD、MDHOVRP 的 扩展 ，MDSVRPPD 的 示意 图 见 图 1(e)。 

分 开 的 研究 ， 将 SVRPPD 与 MDHOVRP 结合 考虑 的 文献 涉 6) 多 中 心 半 开放 式 送 取 需 求 可 拆 分 的 车 辆 路 径 问题 

及 较 少 。 在 现实 复杂 的 物流 配送 网 络 中 ， 仅 靠 解决 以 上 单一 MDHOSVRPPD 问题 中 配送 中 心 具 有 半 开 放 式 的 特点 。 


| 


问题 是 不 全 面 的 ， 而 需要 综合 考虑 车 辆 共享 、 各 配送 中 心 2 开放 式 、 闭 合式 及 半 开 放 式 三 者 是 由 配送 中 心 对 车 辆 始末 停 
间 联 合 配 送 以 及 对 客户 点 送 取 需求 合理 拆 分 的 情况 ， 才 能 更 。 靠 位 置 进行 区 别 的 。 开 放 式 的 配送 中 心 允许 车 辆 在 行驶 路 径 
好 的 优化 现实 物流 配送 路 径 。 中 自由 停靠 ， 车 辆 最 终 不 一 定 要 返回 配送 中 心 ， 闭 合式 的 配 

MDHOSVRPPD 问题 是 现实 物流 配送 模式 发 展演 变 而 来 ，“ 送 中 心 要 求 车 辆 的 始末 位 置 是 配送 中 心 且 必须 一 致 ， 半 开放 
有 多 中 心 、 半 开放 式 以 及 客户 点 送 取 需 求 可 拆 分 等 特点 。 式 的 配送 中 心 要 求 车 辆 完成 任务 后 可 就 近 返 回 任意 配送 中 心 ， 
在 运行 至 客户 点 时 ， 车 辆 负载 减少 或 增加 ， 导 致 负载 波动 ， 车 辆 可 以 共享 配送 网 络 。 因此 MDHOSVRPPD 问题 是 进一步 
必须 对 车 辆 每 个 线路 进行 可 行 性 检查 ， 进 一 步 增加 了 算法 的 降低 配送 成 本 的 演变 ， 也 是 MDSVRPPD 问题 的 拓展 。 
求解 难度 。 为 了 更 好 地 求解 MDHOSVRPPD 问题 , 提升 求解 ” MDHOSVRPPD 的 示意 图 见 图 1(f)。 
效率 和 解 的 质量 , 本 文 设计 了 一 种 针对 MDHOSVRPPD 特点 关于 MDHOSVRPPD 问题 的 最 优 解 理论 ， 文 献 20] 已 证 
的 大 变异 邻 域 遗传 算法 ， 采 用 二 维 染 色 体 编码 ， 不 仅 实现 了 明 多 中 心 半 开放 式 下 的 车 辆 运输 ， 仅 影响 车 辆 在 配送 中 心 的 
对 客户 点 送 取 需 求 的 拆 分 ， 而 且 提高 了 求解 效率 ， 并 引入 顺 首 末 停靠 位 置 ， 而 配送 车 辆 对 客户 点 送 取 货 的 行驶 路 径 不 受 
序 交 叉 算 子 ， 加 强 种 群 之 间 的 交流 则 在 通过 该 算法 以 获得 较 影响。 因此 由 Dror 和 Trudeau 提出 的 需求 拆 分 最 优 解 相关 理 
高 质量 的 满意 解 。 论 ， 也 适用 于 本 文 研究 的 MDHOSVRPPD 问题 。 


1 MDHOSVRPPD 的 提出 


MDHOSVRPPD 问题 是 对 传统 VRP 问题 的 扩展 。 由 以 下 
VRP 问题 扩展 而 来 : 

1) 带 回程 车 辆 路 径 问 题 

带 回 程 车 辆 路 径 问 题 (VRP with Backhauls，VRPB)05 是 
在 配送 车 辆 只 能 进行 后 尾 装载 (Rear-Loaded) 的 状态 下 产生 的 。 
配送 车 辆 负载 货物 从 配送 中 心 出 发 ， 先 服务 具有 送 货 需 求 的 
客户 点 ， 再 服务 具有 取 货 需求 的 客户 点 ， 最 后 车 辆 回 到 配送 
中 心 。VRPB 的 示意 图 见 图 1(a) 。 

2) 混 合 装载 带 回 程 车 辆 路 径 问题 

随 着 物流 运输 技术 的 进步 ,配送 车 辆 不 仅 可 以 后 尾 装载 ， 
还 可 以 进行 侧面 装载 (Side-Loaded)。 因 此 ， 对 送 取 顺 序 不 会 
加 以 限制 ， 产 生 混 合 装 载 带 回程 的 车 辆 路 径 问 题 (VRP with 
Backhauland Mixed Load, VRPBM)Hs1, VRPBM 的 示意 图 见 图 1(b)。 
VRPB 与 VRPBM 问题 都 是 解决 单一 配送 中 心 、 客 户 点 
单一 配送 需求 的 车 辆 路 径 问 题 。 随 着 时 代 的 发 展 解决 现 
有 物流 运输 问题 具有 一 定 的 局 限 性 。 

3) 同 时 送 取 货 车 辆 路 径 问 题 (e)MDSVRPPD (DMDHOSVRPPD 

同时 送 取 货 车 辆 路 径 问 题 (VRP with Simultaneous 全 | as 入 ee 的 EE 


Delivery and Pickup, VRPSDP)0715 是 指 客户 同 时 具有 送 货 、 

取 货 两 种 需求 ， 且 车 辆 对 其 只 能 服务 一 次 ，VRPSDP 问题 是 图 1 MDHOSVRPPD 提出 过 程 

对 VRPBM 问题 的 扩展 。VRPSDP 的 示意 图 见 图 1(c)。 Fig. 1 Presentation process of MDHOSVRPPD 
4) 送 取 需 求 可 拆 分 车 辆 路 径 问 题 


1 凑 尝 模型 
在 实际 应 用 中 ， 部 分 客户 点 需求 量 较 大 ， 如 果 仍 要 求 每 ”2 ”MDHOSVRPPD 的 数学 模型 


个 客户 点 仅 由 一 辆 车 提供 服务 ， 不 仅 造 成 车 辆 空 载 严 重 ， 还 2.1 问题 描述 
兽 加 了 车 辆 使 用 数量 ， 使 派遣 成 本 增加 。Dror 和 Trudeaut5] MDHOSVRPPD 问题 可 以 描述 如 下 :给 定 多 个 配送 中 心 、 
首次 提出 了 送 货 需求 可 拆 分 的 车 辆 路 径 问题 (Split Delivery 客户 点 ， 各 个 配送 中 心 之 间 可 联合 配送 ， 客 户 点 同时 具有 六 
VRP, SDVRP)， 合 理 拆 分 使 得 客户 点 被 配送 车 辆 服务 一 次 的 ” 货 需求 与 取 货 需求 且 单 个 客户 点 的 总 需求 可 能 超过 车 辆 载 
限制 取消 了 。 这 样 不 仅 提高 了 车 辆 的 利用 率 ， 还 可 以 节省 运 同一 类 型 的 车 辆 以 负载 状态 从 配送 中 心 出 发 ， 携 带 客 户 点 
输 距 离 和 服务 的 车 辆 数目 。 后 来 研究 学 者 开始 探讨 对 客户 点 。 需要 的 货物 ， 对 客户 点 执行 先 送 货 、 后 取 货 服务 ， 服 务 时 

送 货 、 取 货 双重 需求 拆 分 的 车 辆 路 径 问 题 (VRP with Split 。 允许 超过 车 辆 最 大 载重 , 车 辆 返回 时 选择 距离 近 的 配送 中 心 。 
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个 客户 点 的 需求 可 以 由 多 辆 车 服务 来 满足 。 目 标 是 找到 


a) 配 送 中 心 、 


上 

组 满足 所 有 客户 点 需求 的 车 辆 路 线 ， 并 使 配送 距离 最 小 化 。 
基于 以 上 分 析 ， 本 文 将 构建 车 辆 运输 距离 最 小 化 的 车 辆 

路 径 模 型 ， 并 满足 以 下 模型 假设 : 

户 点 位 置 已 知 ; 


b) 每 个 客户 点 的 送 取 需 求 量 已 知 ; 
中 所 有 车 辆 起 止 点 均 为 配送 中 心 ， 车 辆 装载 能 力 相 同 且 
运输 过 程 中 实际 装载 量 不 得 超过 车 辆 载重 ; 
的 送 取 总 需求 允许 超过 车 辆 载重 ， 客 户 点 可 
以 接受 多 次 服务 ; 


dd) 各 客户 点 


9) 车 辆 对 客户 点 的 服务 无 顺序 要 求 ; 


人 车辆 对 客 
2.2 ”数学 模型 


户 点 的 服务 无 时 间 要 求 。 


以 距离 最 小 化 为 优化 目标 , 建立 的 MDHOSVRPPD 数学 


模型 如 下 : 


minz=2, 2 2 rnCy (1) 
ieV jeV keK 
> 方 -》 yi=D,;jeLUB (0O) 
ieLUB 
-Sw =D;jeLUB (G3) 
PD 琶 0 (4) 
2 Din- 2 2 rn-0k ek G) 
iEe ieLUB jeLUB 
2 2 r= 2 2 rn = bkek (6) 
ieLUB jeM 
ty QD iiieV (7) 
Xi >0, yy >0,27 >0, x EZ (8) 


目标 函数 式 (1) 为 距离 最 小 化 ， 式 (2)-(3) 保 证 所 有 客户 点 
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炎 
沽 


= 


说 明 


配送 中 心 合集 M e{1,2,L ,m} 
具有 送 货 需求 的 虚拟 客户 点 合集 
Le{fm+l,m+2,L ,m+n} 


具有 取 货 需求 的 虚拟 客户 点 合集 


Be{m+n+t+l,m+n+2,L ,m+2n} 


配送 中 心 和 客户 点 合 旨 
车 辆 的 标号 (车 辆 型 号 统一 ) ke{l,2,L ,K} 
使 用 车 辆 的 总 数 
车 辆 的 最 大 装载 能 
从 配送 中 心 运输 到 客户 点 7 的 总 需求 量 ; 


jeL 


在 客户 点 j 取 货 并 运 回 配送 


jeB 


EV=MULUB 


P 心 的 总 需求 量 ; 


路 径 中 任意 点 i 到 点 j 的 距离 ，i,jeV 
车 辆 从 点 i 行驶 至 点 j 的 次 数 ，i,jeV 
经 过 客户 点 1 ,7 后， 车辆 上 剩余 的 待 送 货物 
量 ; i,jeLUB 


经 过 客户 点 i, j 后 ， 车 辆 上 需 运 回 配送 中 心 的 


货物 量 ;，i,jeLUB 


解码 过 程 中 ， 车 软 
分 为 两 阶段 ， 第 一 阶段 从 左 到 右 将 客 


pop _ Pnum 
取 货 任务 


按照 贪 禁 原 则 ， 车 有 
约束 检验 , 直到 车 辑 


外 足 送 货 企 


j 没 有 待 配送 货物 


的 需求 都 被 满足 ; 式 (4) 确 保 车 辆 离开 配送 中 心 时 车 上 只 有 待 


配送 的 货物 , 当 车 加 


式 (5) 表 示 车 辆 服务 客 


返回 配送 中 心 时 ,车 上 只 有 收取 的 货物 ; 
户 点 的 次 数 等 于 离开 客户 点 的 次 数 ， 车 


辆 不 能 在 客户 点 进 和 
发 且 仅 离开 一 次 ， 
实际 装载 不 超过 车 辆 载 寻 


了 停靠 ; 式 (6) 表 示 车 辆 从 任意 配送 中 心 出 
最 终 返 回 任意 配送 中 心 ， 式 (7) 保 证 车 辆 
EE;， 式 (8) 为 变量 的 取 值 约束 。 


模型 中 各 符号 含义 见 表 1 所 列 。 
3 变异 邻 域 遗 传 算法 


MDHOSVRPPD 是 NP 难 问题 ， 其 送 取 需 求 可 拆 分 的 特 
点 ， 增 加 了 求解 难度 。 对 于 送 取 需 求 可 拆 分 的 车 辆 路 径 问 题 


的 求解 ,大 多 数学 者 采 


构造 型 启发 式 算法 上 62， 本 文 对 于 


MDHOSVRPPD 的 求解 采用 现代 启发 式 算法 ， 现 代 启 发 式 算 


法 对 于 解决 大 规模 的 VRP 类 问题 能 够 在 合理 的 时 间 内 得 出 
近似 最 优 解 。 将 大 变异 遗传 算法 (Genetic Algorithm based on 


big Mutation, GAM) 和 邻 域 算法 Rl(Variable Neighborhood 


Search, VNS) 结 合 


求解 MDHOSVRPPD 问题 ， 前 者 具有 出 色 


的 全 局 搜索 能 力 ， 后 者 有 较 强 的 局 部 搜索 能 力 ， 因 此 本 文 设 
计 了 大 变异 邻 域 遗 传 算 法 (Genetic algorithm based on big 
Mutation and Variable Neighborhood Search, GAMVNS )。 


3.1 编码 与 解码 

本 文采 用 两 段 式 实数 编码 ， 将 真实 客户 点 在 原 位 置 的 者 
础 上 分 成 两 个 虚拟 客户 点 ， 虚 拟 客户 点 只 具有 送 货 任务 或 者 
段 为 送 货 任务 随机 全 排列 ， 第 二 段 为 对 
排列 ， 所 有 任务 储存 在 pop_Pnum 中， 记 


取 货 任务 。 编 码 第 
应 的 取 货 任务 随机 全 1 
录 任 务 的 服务 顺序 如 民 


阶段 根据 每 辆 车 服务 首 
车 的 始末 配送 


E 离 最 短 原 则 ， 确 定 每 各 
到 2 进行 求解 ， 假 设 车 辆 
户 点 3 时 满足 送 货 任务 和 一 部 分 取 货 任务 ， 访 
时 满足 部 分 送 货 任 务 和 全 部 取 货 任务 ， 此 时 车 辆 
日 已 满载 ， 则 退出 路 径 分 配 。 派 出 车 辆 己 
及 务 客户 点 5 的 送 货 任 务 和 客户 点 3 
户 点 1， 假 设 总 满足 客 
则 的 分 配 过 程 和 如 、 厂 的 最 后 路 径 如 图 


没有 竺 配送 货物 ] 


的 取 货 任务 ， 再 去 服务 客户 点 2 和 客 


如 加 加 四国 燥 回 轩 浅 因 
图 2 编码 示意 图 
Fig.2 Coding diagram 
客户 点 服务 执行 先 送 后 装 原则 。 解 码 
户 点 的 任务 划分 给 车 加 
E 务 后 尽 可 能 多 的 取 货 ， 并 进行 


待 运 


习 货 物 满载 时 ， 


分 配 阶段 ， 等 待 返回 配送 中 心 ， 派 出 新 车 对 剩余 


j 服 务 首 


路 径 ， 再 对 车 辆 的 首 


i 


1 
ka 


标 函 数值 。 


标 函 数 式 (1) 确 定 ， 
数 ， 则 染色 体 s 的 适应 度 函 数 上 如下; 


尾 客 户 点 到 配送 中 
Pi 人心 。 


j 满足 客 


户 点 4 的 全 部 


尾 客 户 点 到 配送 中 心 距离 最 小 来 
的 始末 配送 中 心 ， 6、7 为 配送 中 心 。 
在 后 续 操 作 时 只 对 pop_Prnum 信 
果 作 完成 后 依据 配送 车 辆 的 载重 约束 重新 为 每 
配 客户 点 ， 得 到 每 个 车 辆 的 运行 
码 操 作 就 可 以 获得 最 终 路 径 。 只 要 pop_Pnum 
看 由 它 产生 的 最 终 路 径 也 一 定 是 正确 的 ， 
Pnum 的 第 一 阶段 编码 决定 了 后 续 编 码 。 通 过 此 编码 方 
点 需求 的 拆 分 ， 而 且 使 得 后 续 扰 动 操 作 
可行 解 ， 提 高 了 算法 效率 。 


息 的 第 一 段 进 


于 构建 距离 最 


(9) 


录用 定稿 


3.3 选择 交叉 操作 


本 文选 择 操作 采用 
方式 。 轮 盘 赌 是 个 体 以 一 定 概率 被 选中 ， 


张 疾 钰 ， 等 : 多 中 心 半 开 放 式 送 取 需求 可 拆 分 的 车 辆 路 径 研究 


选中 的 概率 则 大 , 反之 则 小 。 精 英 保留 策略 


体 的 适应 度 函 数值 大 小 确定 ， 个 体 适 应 


CSET 
2 回回 四 加 
52 回国 四 
2222 回 四 加 


轮 盘 赌 与 精英 保留 策略 相 结合 的 


概率 大 小 根据 个 
度 函 数值 高 的 被 


pop _ Pnum 


光 国 回国 四 
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应 度 函 数 


第 
是 适 


大 | 


个 


体 间 基因 充 


1 :分配 过 程 4p 一 433p 一 5p 一 5p 


首尾 解码 


| 


ki:i6—p4,— 4 3 3 5 7 


>6 


3 2 
2 slalil < 
315121 1 
ZA 3 VA? | 1 
图 3 解码 示意 
Fig.3 Decoding diagram 


图 


] :7 一 上 7 


p32 2p 


>1,—>1, 


将 个 体 中 的 客户 点 3 与 客 


值 最 高 的 父 代 代替 适应 度 函 数值 最 低 的 子 代 。 采 取 这 两 
种 策略 有 助 于 加 强 所 有 
成 稳定 的 下 一 代 。 


分 交流 以 及 快速 形 


个 邻 域 操作 。 直 至 循环 到 最 后 一 个 
止 输出 最 优 解 。 
5(a)， 将 个 体 中 的 客 


户 点 3 移动 到 客户 


户 点 6 


交叉 操作 则 是 采用 顺序 交叉 算 子 。 对 个 体 pop_1 顺序 交叉 。 达到 WMax_se 时 ， 进 入 下 
时 ， 从 种 群 中 随机 选取 Pop -2 个体 ， 分 别 随机 选取 pop_1、 邻 域 操 作 ， 搜 索 终 
pop_2 中 四 个 基因 gs 、sg 、sg 和 8 ，pop_1 中 的 sg、s 之 a) 插 入 规则 如 图 
间 的 基因 作为 newpop_1 个 体 的 第 一 段 , 消除 pop -2 中 su 、so 之 点 6 后 产生 新 的 个 体 。 
间 的 客户 点 ， 保 持 剩 下 客户 点 位 置 不 变 ， 作 为 newpop_1 的 第 二 b) 交 换 规 则 如 图 5(b)， 
段 , 此 时 newpop_1 个 体 生产 完毕 , newpop_2 同 理 。 如 图 4 所 示 。 进行 交换 产生 新 的 个 体 。 

的 人 c)2-opt 规则 如 图 
wa 间 的 客户 逆序 排列 ， 客 
序 排列 ， 从 而 产生 新 的 个 体 。 


4 


顺序 交叉 示意 


Fig.4 Sequence crossing diagram 


3.4 大 变异 操作 


传统 遗传 算法 容易 “早熟 ” 为 了 让 种 群 
陷 。 大 变异 策略 具 
， 并 记录 每 一 代 种 群 的 最 高 适应 


采取 大 变异 策略 弥 # 


设 定 一 个 合适 的 变异 率 搜索 


度 所. 和 平均 适应 度 


这 一 缺 


J > 当 满足 : 


OX frr < fa 


表示 该 代 种 群 


异 率 大 5 倍 以 上 的 概率 对 该 代 种 群 进 


0.5<a<1，0 为 密集 


opt, 来 增强 算法 的 


和 最 大 搜索 次 数 Max_ 


人 


大 


四 主 h 2 
4 早熟 a 


说 


A 


有 高 度 “ 集 中 ”的 表现 ， 随 后 以 普通 


体 为 P31， 首先 


(10) 


子 ， 表 示 种 群 的 集 
变异 操作 采用 了 三 种 邻 域 结 构 ， 


行 


次 变异 操作 。 其 
程度 。 


分 别 是 插入 、 交 换 、2- 


5(c)， 将 个 体 中 客户 点 3 与 客户 点 6 之 
户 点 3 保持 不 变 ， 客 户 点 5、2、6 道 


(b) 交换 


Ass oA |? 


受 | 


5 


邻 域 搜索 示意 


(c) 2-opt 


多 


Fig.S$ Neighborhood search diagram 


Se 次 最 优 解 不 变 ， 进 


4 


原 概率 变异 


实验 与 结果 分 析 
本 文 依次 选取 MDVRP 算 例 、MDHOSVRPPD 算 例 与 


新 一 代 种 群 


N 


局 部 搜索 。 设置 邻 域 结构 的 最 优 解 次 数 se 使 用 三 种 策略 7 
Se ， 当 个 体 进行 第 一 个 邻 域 操作 时 , 在 
入 下 一 个 邻 域 操作 ， 或 当 个 体循环 次 数 。” 率 。 算 法 流程 图 如 图 6。 


2-opt 
N 


6 算法 流程 
Fig.6 Algorithm flow 


SVRPPD 算 例 验证 本 文 算法 
仅 证 明了 本 文 算法 的 有 


chart 


邻 域 解 改进 


生 新 个 体 ， 能 大 大 提高 个 体 的 适应 
进而 提升 GAMVNS 算法 的 局 部 搜索 能 力 ， 提 高 算法 运 


度 ， 
算 效 


的 有 效 性 。3 


其 中 SVRPPD 算 例 不 


效 性 


， 而 且 证 明 


半 


-放送 取 和 需求 可 于 


Ff 分 的 配送 模式 优 于 


了 联合 配送 下 多 中 心 
单 中 心 下 送 取 需 求 可 


| ee 
ChninaXIv 
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拆 分 的 配送 模式 。 为 了 测试 GAMVNS 算法 的 性 能 ， 对 算 例 表 5 实验 二 算 例 数据 
进行 大 量 实验 验证 。 本 文 在 MATLAB2016B 中 编程 ， 并 使 用 Tab.5 Experiment2 example data 
Inter®@CoreT™Mi5-1035G1 CPU 2.19GHz 计算 机 系统 在 Windows 编号 xi yi di pi 
10 上 执行 ， 该 系统 具有 16GB 内 存 。 实 验 基本 参数 设置 见 表 1 37 52 5 2 
2，n、N、G 、Pc、Pm、Se 、Max_Se 分 别 为 客户 规模 、 种 2 49 49 30 0 
群 规模 、 最 大 迭代 次 数 、 交 叉 概率 、 变异 概 率 、 最 优 解 次 数 、 3 52 64 13 3 
最 大 邻 域 搜索 次 数 。 4 20 26 7 2 
表 2 参数 设置 5 40 30 16 5 
Tab. 2 Parameter setting 6 21 47 7 8 
n N G Pc Pm Se Max_ Se 7 17 63 S 14 
0<n<50 100 50~500 0.8 0.08 50~100 70~100 8 31 62 12 12 
S50<nz75 300 500~700 0.8 0.08 50~100 70~100 9 52 33 7 4 
75<m<100 500 700~1000 0.8 0.07 50~100 70~100 10 S51 21 
4.1 MDVRP 算 例 验证 11 42 41 19 0 
实验 一 为 了 验证 本 文 算法 的 有 效 性 ， 选 取 文 献 [23~25] 12 31 32 28 1 
中 的 MDVRP 算 例 进 行 验证 。 该 算 例 中 包括 30 个 客户 点 ( 编 13 5 25 5 18 
号 为 1-30), 3 个 配送 中 心 (编号 为 31-33), 配送 车 辆 的 最 大 容 14 12 42 15 
量 为 10t 及 最 大 行驶 距离 为 50km， 所 有 客户 点 需求 都 在 2t 15 36 16 6 
及 以 下 。 表 3 给 出 文献 [23] 的 狼 群 算法 (WPA)、 文 献 [24] 的 竞 16 52 41 12 3 
争 决策 算法 (CDA) 以 及 文献 [25] 的 禁忌 搜索 算法 (TS)、 量子 遗 17 27 23 3 0 
传 算法 (QGA) 与 本 文 提 出 的 GAMVNS 算法 运行 10 次 的 结果 18 17 33 21 20 
对 比 。 表 中 KOS 表示 已 知 最 优 解 ，Best 、Worst 、Avg 分 别 表 19 13 13 9 0 
示 算 法 运行 10 次 得 到 的 最 优 解 、 最 差 解 以 及 平均 解 , %Dev 、 20 57 58 28 0 
CPU 为 最 优 解 偏 差 与 算法 平均 运行 时 间 。 配送 车 辆 行驶 路 径 21 62 42 5 3 
见 表 4。 22 42 57 6 2 
表 3 可 以 看 出 ,本文 算法 求 得 最 优 解 与 已 知 最 优 解 23 16 57 4 12 
相同 均 是 116.01， 而 求 得 的 最 差 解 及 平均 解 均 优 于 其 他 24 8 52 2 8 
算法 所 求 。 从 求解 时 间 方 面 来 看 ， 本 文 算法 运行 10 次 的 25 7 38 5 23 
平均 时 间 为 7.84s， 运 行 时 间 比 其 他 算法 有 显著 降低 ， 有 26 27 68 3 4 
较 强 的 优势 。 27 30 48 9 6 
表 3 实验 一 结果 对 比 28 43 67 9 5 
Tab.3 Comparison of experiment 1 results 29 58 48 5 1 
算 例 ”KOS 总 距离 总 距离 总 距离 pe 30 58 27 9 10 
Best %Dev Worst %Dev Avg %Dev 31 37 69 6 5 
WPADP3] 122.42 5.53 - - - - - 32 38 46 10 2 
CDAR4 123.33 6.31 - - - - - 33 46 10 5 18 
TSP5 116.01 116.01 0 343.31 195.93 187.34 61.49 243 34 61 33 14 12 
QGAP5] 116.01 0 326.48 181.42 176.37 52.03 216 35 62 63 17 0 
GAMVNS 116.01 0 150.13 29.41 142.81 23.10 7.84 36 63 69 1 
表 4 实验 一 最 优 解 路 径 37 32 22 6 3 
Tab.4 Experiment 1 optimal solution path 38 45 35 12 3 
算法 车 辆 行驶 路 线 总 距离 39 59 15 10 
33-6-3-18-26-12-5-25-33 40 5 6 6 1 
33-17-21-19-20-23-24-2-9-27-16-33 41 10 17 16 11 
GAMVNS 116.01 
32-11-29-18-13-8-15-1-32 42 21 10 6 7 
31-4-7-10-22-14-30-31 43 S 64 1 10 
4.2 MDHOSVRPPD 算 例 验证 44 30 15 8 8 
实验 二 MDHOSVRPPD 问题 的 相关 约束 较 多 ,有 多 个 配 45 39 10 3 7 
送 中 心 、 多 个 客户 点 、 客 户 点 送 取 需 求 等 特点 。 目 前 没 46 32 39 4 1 
通用 的 算 例 集 ， 本 文选 用 标准 算 例 库 中 的 MDVRP 算 例 47 25 32 20 5 
P02( 算 例 集 来 源 https://neo.lcc.uma.es/vrp/vrpinstances/multiple- 48 25 55 8 9 
depot-vrp-instances/), 并 按 MDHOSVRPPD 问题 的 特点 修改 后 49 48 28 11 8 
来 进行 验证 。P02 算 例 包括 50 个 客户 点 (编号 为 1-530)，4 个 50 56 37 7 3 
配送 中 心 (编号 为 51-54)， 和 车辆 容量 为 160。 客 户 的 送 取 需求 51 20 20 0 0 
由 NagyEg9 提 出 的 需求 分 离 规则 产生 : 坟 、 六 为 客户 点 的 化 52 30 40 0 0 
标 ，d; 、pi 为 客户 点 送 货 、 取 货 需 求 ，g; 为 客户 点 原始 需求 ， 53 50 30 0 0 
5=min 712 1) 计算 每 个 客户 比率 ， 由 d=gn 和 p=9.(1-i) 54 60 50 0 0 
计算 得 出 各 个 客户 点 的 送 取 需 求 ， 算 例 具 体 信息 见 表 5。 计算 结果 稳定 。 并 与 传统 遗传 算法 (Genetic Algorithm， 
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GA)、 传 统 蚁 群 算法 (Ant Colony Optimization，ACO)、 传 统 4.3 配送 模式 对 比 验证 
大 变异 遗传 算法 (Genetic Algorithm based on big Mutation, 实验 三 为 验证 联合 配送 模式 下 多 中 心 半 开放 式 送 取 需 
GAM)、 传 统 变 邻 域 搜索 算法 (Variable Neighborhood Search， ” 求 可 拆 分 配送 模式 优 于 独立 配送 模式 下 单 中 心 送 取 需 求 可 拆 
VNS) 作 对 比 ， 结 果 见 表 6 和 表 7，Best 、Num 为 算法 运行 20 分 配送 模式 (VRP with Split Deliveries and Pickups, SVRPPD)， 
次 的 最 优 解 和 车 辆 数 ，%Dev 、cPU 为 最 优 解 偏差 与 算法 平均 。 选取 了 MDVRP 标准 算 例 库 中 的 P01、P02、P03、P04、P05、 
运行 时 间 。 P06 进行 半 开 放 式 送 取 需 求 可 拆 分 的 多 中 心 与 单 中 心 对 比 看 
由 表 6 看 出 ， 传 统 遗 传 算法 、 大 变异 遗传 算法 、 蚁 群 算 ” 究 。 多 中 心 的 配送 中 心 坐 标 取 算 例 中 的 数据 ， 单 中 心 的 配送 
法 及 变 邻 域 搜 索 算 法 搜索 策略 单一 , 容易 出 现 * 早熟” 现象 ， 中 心 坐 标 见 表 8， 客 户 点 送 取 需求 数据 的 处 理 按照 实验 二 的 
在 求解 质量 上 表现 较 差 。 本 文 提 出 的 GAMVNS 算法 与 四 种 。 标准 执行 。 表 9 中 Best 、Worst 、Ahvg 分 别 为 算法 运行 20 次 
算法 相 比 ， 在 配送 总 距离 方面 分 别 优 化 了 5.66%、6.88%、 得 到 的 最 优 解 、 最 差 解 以 及 平均 解 。 %Dev 为 最 优 解 偏差 ， 
2.83%、16.46% ， 求 解 质量 方面 均 优 于 传统 算法 。 从 求解 时 Num 为 车 辆 数 ， CPV 为 算法 运行 20 次 的 平均 时 间 。 结 果 见 
间 上 看 ， 本 文 算法 的 平均 时 间 为 7.2s， 计 算 时 间 比 较 理想 。 表 8 和 表 9。 
最 优 解 偏差 为 0.8% ,能够 有 效 的 求解 出 最 优 解 。 表 7 为 配送 表 9 看 出 ， 联 合 配送 模式 与 独立 配送 模式 下 ， 各 算 例 
车 辆 行驶 路 径 。 由 图 7 的 最 优 路 径 的 迭代 图 可 以 看 出 ， 本 文 ”求解 结果 中 , 最 优 解 偏差 最 大 为 8.9%, 最 小 为 0.5%, 求解 时 


Pal 


— 


的 GAMVNS 算法 能 够 以 较 少 的 迭代 次 数 快速 地 进行 收敛 ， 间 均 在 8s 内 。 可 以 看 出 ,本文 算法 的 优化 结果 求解 速度 较 快 ， 
有 较 强 的 寻 优 能 力 。 可 以 稳定 的 收敛 到 较 优 解 。 
表 6 实验 二 结果 对 比 对 比 联合 配送 模式 下 多 中 心 半 开 放 式 送 取 需求 可 拆 分 的 
Tab.6 Comparison of experiment 2 results 配送 模式 与 独立 配送 模式 下 单 中 心 送 取 需求 可 拆 分 的 配送 模 
算法 Best Num %Dev CPU/s 式 , 可 以 看 出 , 配送 距离 对 应 缩短 为 31.86%、13.0%、25.57%、 
GAMVNS 738.59 3 0.83 7.2 9.89%、27.95% 和 6.55%。 从 表 9 与 图 8 可 以 得 出 ， 联 合 配送 下 
GA 833.41 3 1.47 11.7 多 中 心 半 开 放 式 送 取 需求 可 拆 分 的 配送 模式 ， 相 对 于 独立 配送 
ACO 860.21 3 1.52 12.1 下 单 中 心 送 取 需求 可 拆 分 模式 ， 前 者 大 幅度 减少 了 车 辆 频繁 往 
GAM 780.44 3 0.87 8.3 返 配送 中 心 产 生 的 配送 路 径 ， 进 而 降低 车 辆 的 配送 成 本 。 
VNS 789.45 3 0.91 11.3 表 8 实验 三 算 例 数据 
表 7 实验 二 最 优 解 路 径 Tab.8 Experiment three example data 
Tab.7 Experiment2 optimal solution path 算 例 数 据 MDHOSVRPPD SVRPPD 
算法 车 辆 行驶 路 线 总 距离 算 例 ”客户 数 ” 车 辆 载重 ”配送 中 心 数 ”配送 中 心 数 配送 中 心 坐标 
54-36-3-28-31-8-26-7-23-24-14-25-18-4-47-12- P01 50 80 4 1 (20,20) 
46-27-1-32-22-48-6-14-25-13-51 P02 50 160 4 1 (20,20) 
GAMVNS 54-29-21-50-9-49-5-38-16-34-30-10-39-33-45- 738.59 P03 3 140 5 1 (40,40) 
15-44-37-42-41-13-40-24-43-52 P04 100 100 2 1 (35,20) 
54-35-20-2-11-32-1-22-48-6-43-17-19-40-51 P05 100 200 2 1 (15,35) 
P06 100 100 3 1 (15,20) 
1500 
i 5 ”结束 语 
1300 上 本 文 针对 联合 配送 下 多 中 心 半 开放 式 送 取 需 求 可 拆 分 的 
车 辆 路 径 问题 进行 了 以 下 研究 。(1) 构 建 了 总 配送 距离 最 短 的 
长 多 中 心 半 开 放 式 送 取 需求 可 拆 分 的 车 辆 路 径 数 学 模型 。(2) 设 
度 2000 计 了 大 变异 邻 域 遗 传 算法 对 其 求解 ， 大 变异 策略 和 邻 域 搜索 
900 十 | 通过 三 
组 算 例 对 模型 和 算法 进行 了 验证 分 析 。 结 果 表 明 ， 本 文 建立 
的 模型 和 算法 能 够 有 效 地 解决 MDHOSVRPPD 问题 , 且 搜 索 
a 稳定 ， 具 有 较为 明显 的 优势 。 
0 100 200 300 400 500 600 700 800 900 a i held le 
迭代 次 数 业 提 供 更 贴 合 现实 的 配送 路 径 ,不 仅 要 考虑 配送 距离 最 小 ， 
还 要 考虑 客户 满意 度 问 题 。 因 此 ， 在 时 间 窗 约束 下 求解 多 中 
图 7 最 短路 径 迭 代 图 心 半 开 放 式 送 取 需 求 可 拆 分 的 车 辆 路 径 问 题 是 下 一 步 的 研究 
Fig.7 Shortest path iteration graph 方向 。 
表 9 实验 三 结果 对 比 
Tab.9 Comparison of the results of the third experiment 
MDHOSVRPPD SVRPPD 
Best Worst Avg %Dev Num CPU/s Best Worst Avg %Dev Num CPU/s 
959.62 1122.08 987.99 2.9 6 7.2 1265.32 1573.06 1293.6 2.2 6 7.3 
738.59 833.41 779.86 5.3 3 7.2 834.60 992.23 864.48 3.5 3 6.9 
917.82 1143.7 1008.58 8.9 6 6.8 1152.52 1416.46 1208.07 4.6 6 7.1 
1955.21 2016.45 1989.43 1:7 9 7.5 2148.69 2241.41 2183.12 1.6 9 | 
1522.29 1642.5 1601.31 4.9 S$ 7:3 1947.81 2057.52 1974.55 1.4 S 6.8 
1804.34 1837.48 1875.6 3.8 9 6.9 1922.45 1924.80 1923.63 0.5 9 7.2 
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图 8 MDHOSVRPPD 配送 模式 下 最 优 路 径 图 


Fig.8 Optimal path diagram under MDHOSVRPPD distribution mode 
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